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Abstract. The origin of galactic and extra-galactic magnetic fields is an unsolved problem in modern 
cosmology. A possible scenario comes from the idea of these fields emerged from a small field, a seed, 
which was produced in the early universe (phase transitions, inflation, ...) and it evolves in time. 
Cosmological perturbation theory offers a natural way to study the evolution of primordial magnetic 
fields. The dynamics for this field in the cosmological context is described by a cosmic dynamo like 
equation, through the dynamo term. In this paper we get the perturbed Maxwell's equations and 
compute the energy momentum tensor to second order in perturbation theory in terms of gauge invari- 
ant quantities. Two possible scenarios are discussed, first we consider a FLRW background without 
magnetic field and we study the perturbation theory introducing the magnetic field as a perturbation. 
The second scenario, we consider a magnetized FLRW and build up the perturbation theory from this 
background. We compare the cosmological dynamo like equation in both scenarios. 
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1 Introduction 

Magnetic fields have been observed on several scales in the universe. Galaxies and clusters of galaxies 
contain magnetic fields with strengths of ~ 10~ 6 G [1, 2], fields within clusters are also likely to 
exist, with strengths of comparable amplitud [3, 4]. There is also, some evidence of their presence 
on scales of superclusters [5-8]. On the other hand, the possibility of cosmological magnetic field has 
been addressed comparing the CMB quadrupole with one induced by a constant magnetic field (in 
coherence scales of ~ IMpc), constraining the field amplitude to B < 6.8 x 10 _9 (ri m /i 2 ) 1 / 2 Gauss [9- 
11] . However, the origin of such large scale magnetic fields is still unknown. These fields are assumed to 
be increased and maintained by dynamo mechanism, but it needs a seed for the mechanism takes place 
[12]. Astrophysical mechanisms, as the Biermann battery have been used to explain how the magnetic 
field is mantained in objects as galaxies, stars and supernova remnants [13-16], but they are not likely 
correlated beyond galactic sizes [17]. It makes difficult to use only such mechanism to explain the origin 
of magnetic fields on cosmological scales. Possible explanations appear in diffcrents context to generate 
the seed from the which the astrophysical mechanism start. For example, magnetic fields could be ge- 
nerated during primordial phase transitions (such as QCD, the electroweak or GUT), parity- violating 
processes which generates magnetic helicity or during inflation [18-27]. Magnetic fields also are gene- 
rated during the radiation era in regions that have nonvanishing vorticity, this seed was proposed by 
Harrison [28-30]. Magnetic fields generation from density fluctuations in prcrccombination era has 
been investigated in [31]. The advantage of these primordial processes is that they offer a wide range 
of coherence lengths (many of which are strongly constrained by Nucleosynthesis [32-34]), while the 
astrophysical mechanisms produce fields with characteristic scales of the size of astrophysical object. 
Recently has been found a lower limit of the large scale correlated magnetic field, this fact constrains 
models for the origin of cosmic magnetic fields, giving a possible evidence for their primordial origin 
[35-37]. 

Cosmological perturbation theory [38-40] is a powerful tool to understand the present properties of 
the large-scale structure of the Universe and their origin. It has been mainly used to predict effects on 
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the temperature distribution in the Cosmic Microwave Background (CMB) [41, 42] furthermore, linear 
perturbation theory combined with inflation suggest us that primordial fluctuations of the universe are 
adiabatic and Gaussian [43]. However, due to the high precision measurements reached in cosmology, 
higher order cosmological perturbation theory is required to test the current cosmological framework 
[44]. There are mainly two approaches to study higher order pcrturbative effects: one uses nonlinear 
theory and different manifestations of the separate universe approximation, using the AN formalism 
[45, 46], and the other one is the Bardeen's approach where metric and matter fields are expanded 
in a power series [47]. Within Bardeen's approach, a set of variables are determined in such way 
that have no gauge dependence. These are known in the literature as gauge-invariant variables 
which have been widely used in different cosmological scenarios [48, 49]. One important result of 
cosmological perturbation theory is the coupling between gravity and electromagnetic fields, finding a 
magneto-geometrical interaction that can change the evolution of the fields on large scales. An effect 
is the amplification of cosmic fields. Indeed, large scale magnetic fields in perturbed spatially open 
FLRW models decay as a -1 , a rate considerably slower than the standard a -2 [50-53], the hyperbolic 
geometry of these open FLRW models lead to the superadiabatic amplification on large scales [54]. 
The main goal in this paper is to study the late evolution of magnetic fields which were generated 
in early stages of the universe. We use the cosmological perturbation theory following the Bardeen's 
approach to find the perturbed Maxwell equations up to second order and also we get an equation 
like dynamo written in terms of gauge invariant variables to first and second order. Futhcrmore, we 
discussed the importance that the curvature and the gravitational potential play in the evolution of 
these fields. The paper is organized as follows: in the next section we briefly give an introduction of 
cosmological perturbation and we address the problem gauge in this theory. In section 3, we write 
the matter equations in the homogeneous and isotropic universe, which was used to generate the first 
and second order dynamical equations. In section 4, we define the first order gauge invariant variables 
for the perturbations not only in the matter (energy density, presure, magnetic and electric field) 
but also in the geometrical quantities (gravitational potential, curvature, shear ..). The first-order 
perturbation of the Maxwell's equations is reviewed in section 5 and together with the the Ohm's 
law allow us to find the cosmological dynamo equation to describe the evolution of the magnetic 
field. The derivation of second-order Maxwell's equations is given in section 7, and following the 
same methodology for the first-order case, we find the cosmological dynamo equation at second order 
written in terms of gauge invariant variables. In the section 8, we use an alternative approximation to 
the model considering a magnetic field in the FLRW background, we find that amplification effects of 
magnetic field appear at first order in the equations, besides of the absence of fractional orders. Also 
we discuss the differences found between both approaches. The final section 9 is devoted to discuss, 
the main results and the connection with future works. 



2 The gauge problem in perturbation theory 

Cosmological perturbation theory help us to find approximate solutions of the Einstein field equations 
through small desviations from an exact solution [55]. In this theory one works with two differents 
space-times, one is the real space-time (M,g a p) which describes the perturbed universe and the other 
is the background space-time (.Mo, g^l) which is an idealization and is taken as reference to generate 
the real space-time. Then, the perturbation of any quantity r (e.g., energy density /x(x, t), 4- velocity 
u a (x,t), magnetic field B l (x,t) or metric tensor g a p) is the difference between the value that the 
quantity r takes in the real space-time and the value in the background at a given point. 1 Now to 
determinate the perturbation in this quantity P, we must have a way to compare r (tensor on the real 
space-time) with (being the value on A4o), this requires an assumption to identify points of 
M. with those of A4o- This is accomplished by assigning a mapping between these space-times called 
gauge choice given by a function X : A4o(p) — > A4(p) for any point p G Aio and p € Ai, which 
generate a pull-back 

X* : M _^ Mo , 9n 
T*(JS) T*(p) ' y > 

1 This difference should be taken in the same physical point. 
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thus, points on the real and background space-time can be compared through of X. Then, the 
perturbation for F is define as 

6r(p) = r(p)-rM(p), (2.2) 

We see that the perturbation ST is completely dependent of the gauge choice because the mapping 
determines the representation on Mq of r(p). However, one can also choice another correspondence y 
between these space-times so that y : .Mo(g) — > -M(p), (p ^ q). 2 A change of this identification map in 
the literature is called gauge transformation. The freedom to choose between different correspondences 
is due to the general covariance in General Relativity, which say that there is no preferred coordinate 
system in nature [56, 57]. Hence, this freedom will generate an arbitrariness in the value of ST at any 
space-time point p, it is called gauge problem in the general relativistic perturbation theory and has 
been treated by [47, 58, 59]. This problem generates unphysical degree of freedom to the solutions in 
the theory, therefore one should fix the gauge or build up quantities that they do not dependent of 
this degree of freedom. 

2.1 Gauge transformations and gauge invariant variables 

To define the perturbation to a given order, we need to introduce the concept of Taylor expansion on 
a manifold and thus the metric and matter fields are expanded in a power series. Following [60-62], 
we consider a family of four-dimensional submanifolds Ai\ with A e M, embedded in a 5-dimcnsional 
manifold N = M x R. Each submanifold in the family represents a perturbed space-time and the 
background space-time is represented by the manifold Mq (A = 0), and on these manifolds we consider 
that the Einstein field and Maxwell's equations are satisfied 



E[ ffA ,T A ]=0 and M[F A ,J A ]=0; 



(2.3) 



each tensor field A, on a given manifold A4\ is extended to all manifold AT through r(p, A) = -T A (p) to 
any p 6 A4\ likewise the above equations are extended to A/". 3 We need to use a diffcomorphism such 
that the difference in the right side of eq. (2.2) can be done. We introduce an one-parameter group 
of diffeomorphisms X\ which identifies points in the background with points in the real space-time 
labelled with the value A. Each X\ is a member of a flow X on M and it specifies a vector field X with 
the property X 4 = 1 everywhere (transverse to the A4 A ) besides other properties, 4 then points which 
lie on the same integral curve of X is to be regarded as the same point [59]. Therefore, according to 
the above, we get a definition for the tensor perturbation 



A A 



r . 



(2.4) 



At higher orders the Taylor expansion is given by [61], 

k=0 

where 



d k x*r 



Now, rewriting eq. (2.4) in the following form 



V— s (k) r 
^ k\ x ' 

k=l 



X=0,Mo 



(2.5) 



(2.6) 



r + X5 { ^r + -5^r 



o(A 3 ), 



(2.7) 



2 This is the active approach where transformations of the perturbed quantities are evaluated at the same coordinate 
point. 

3 In eq. (2.3), g\ and T A are the metric and the matter fields on M\, similarly F\ and J\ are the electromagnetic 
field and the four-current on M\. 

4 Here we introduce a coordinate system x a through a chart on M\ with a = 0, 1, 2, 3, thus, given a vector field on 
J\f, it's the property that X 4 = A in this chart, while theirs other components remains arbitrary. 
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we see in the eqs. (2.6) and (2.7) that representation of r on Aio is splitting in the background value 
To plus 0(k) perturbations in the gauge X\. Therefore, the fc-th order O(k) in r depends on gauge 
X. With this description the perturbations are fields living in the background. The first term in eq. 
(2.4) admits an expansion around of A = given by [61] 

*tr\ Mo = -^c k x r\ Mo = ex P (\c x )r\ Mo , (2.8) 

fe=0 

where CxT is the Lie derivative of r with respect to a vector field X which generates the flow X. 
If we define X?r\jvi = and doing the same for other gauge choice y, using eqs. (2.4)-(2.8), the 
tensor fields F x can be written as 



\ k \ k 

^ = E^4 fc) ^ = E^iM„, (2.9) 

k=0 k=0 



00 \ k 00 \k 

^ = E^4^ = E^i-o, (2.10) 

fe=0 ' k=0 



if r* = r? for any arbitrary gauge X and y, we say that r is totally gauge invariant. We can also 



say that r is gauge invariant to order n ^ 1 iff satisfy 5y T = S^T, or in other way 



A A 

y 1 - u x 

£ x 6 {k) r = 0, (2.11) 



for any vector field X and Vfc < n. To first order (k = 1) any scalar that is constant in the background 
or any tensor that vanished in the background are gauge invariant. This result is referred to as the 
Stew art- Walker Lemma [63], i.e., eq. (2.11) generalizes this Lemma. However, when r is not gauge 
invariant and we have two gauge choices X\ and y\, the representation of r\M i s different depending 
of the gauge used. To transform the representation from a gauge choice X x *r\M a to another y x r\M a 
we use a map <&\ : Mo —> Mo given by 

$ A ee x_ x o y x = $* x r? , (2.12) 

thereby, the diffeomorphism $> x induce a pull-back $^ which changes the representation of r 
in a gauge X\ to the representation of r in a gauge y x . Now, following [64] and using the 
Baker-Campbell-Haussdorf formula [65], we can generalize eq. (2.8) to write ^x^x m the following 
way 

** x r? = expf^ ¥ £ a Jrf, (2.13) 
where is any vector field on Mx- Substitute eq.(2.13) in eq.(2.12), we have explicitly that 

ry = r* + M^r* + y (4 + c i2 ) r* + o(\ 3 ). (2.14) 

Thus, if substitute eq.(2.9) and eq.(2.10) into eq.(2.14), the relation to first and second order pertur- 
bations of r in two different gauge choice are given by 

5$ ) r-8% ) r = c tl r , (2.15) 
s^r - 5fr = 2C (l s^r + (4 + c i2 ) r , (2.16) 

where the generators of the gauge transformation $ are 

t!=Y-X and G2 = [X,Y]. (2.17) 

The gauge transformation given by the eqs. (2.15) and (2.16) are quite general, to first order we 
see that r is gauge invariant iff A^-fo = 0, while to second order one must have another condition 
C^5^ r = 0, and so at high orders. We will apply the formalism described above to the Robertson- 
Walker metric, where k does mention to the expansion order. 
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3 FLRW background 



We assume that at zero order (background) , the universe is well described by a spatially flat Fricdman- 
Lemaitre-Robertson- Walker metric (FLRW) 

ds 2 = a 2 (r) {-dr 2 + 6, ,<!.,■■ d.r' ) , (3.1) 

with a(r) the scale factor. Hereafter the Greek indices run from to 3, and the Latin ones run from 1 
to 3, we will work with conformal time r, and a prime denotes the derivative with respect to r. The 
Einstein's components tensor in this background are given by 

Gq = -jp, (3.2) 

3 = -i( 2 £ (3.3) 
G° = G l = 0, (3.4) 

with H = 2- the Hubble parameter. We consider the background filled with a single barotropic fluid 
where the energy momentum tensor is 

T UD „ = (M(0) + P ( o)) «? 0) «<°> + P(0)K, (3-5) 

with /i( ) the energy density and P/ ) the pressure. The comoving observers are defined by the 
four-velocity u v = (a -1 , 0, 0, 0) with u u u v = —1 and the conservation law for the fluid is 

// (0) +3ifGu (0 )+P(o)) = 0. (3.6) 

To deal with the magnetic field, the space-time under study is the fluid permeated by a weak magnetic 
field, 5 it is a stochastic field and can be treat as a perturbation on the background [66, 67]. Due to 
the magnetic field has no background contribution, the electromagnetic energy momentum tensor is 
automatically gauge invariant (see eq.(2.15)). Now we use the spatial part of Ohm's law [68], which 
is the current proyected in the slices 

(9fii + u^Ui) f = ag\ i g ai _ l F Xa u' 1 , (3.7) 

where j M = (g, J 1 ) is the 4-current and F Xa is the electromagnetic tensor given by 



F = 



a? (r) 






E' 


EJ 


E k \ 


E i 





B k 


—Bi 


E j 


_ B k 





B l 


E k 


Bi 


-B l 


o 1 



(3.8) 



At zero orden in eq. (3.7) we find the usual Ohm's law which gives us the relation between the 
3-current and the electric field 

Ji=<rEi, (3.9) 

where a is the conductivity. We work under MHD approximation, thus in large scales the plasma is 
globally neutral and charge density is neglected (g = 0) [3]. If the conductivity is infinite (er — > oo) 
in the early universe [69, 70], then eq. (3.7) say us that electric field must vanish (Ei = 0) in order 
to keep the current density finite [71, 72]. Hovewer, the current also should be zero (Ji = 0) because 
a nonzero current involves a movement of charge particles which breaks down the isotropy in the 
background. 



3 With the property B? Q s < /f( ) [66]. 
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4 Gauge invariant variables at first order 



We write down the perturbations on a spatially flat Robertson- Walker background. We shall consider 
the fluctuations of the metric and matter fields. The pcrturbative expansion at k— th order of the 
matter quantities arc given by 



oo ^ 

M = M(o) +X!lT /i ( fe )' 
fc=l 

oo 1 

fc=i 

^(g^))> 



(4-1) 
(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 
(4.8) 



Where the fields used in above formulae are the average ones (i.e. B 2 = (B 2 )) 6 We also consider 
the perturbations about a FLRW background, so that the metric tensor is given by 



goo = -a( 



^) 2 ( 1 + 2 E^ (fe) ) 

oo _^ 

9oi = a{r) 2 -n 



k=l 



.,(k) 
hi"* ' 



m = «( r ) 2 



oo \ oo 1 



(fc) 



fc=i 



fe=i 



(4.9) 
(4.10) 
(4.11) 



The perturbations are splitting into scalar, transverse vector part, and transverse trace-free tensor 

;f)=^ fc )ll+ W f^, (4.12) 



with 9 l w- fe ^ = 0. Similarly we can split as 



(fc) n fife) II , a ( fc )-L , o ( fe )-L , (fe)T 



(4.13) 



3 This happens because the average evolves exactly like B 2 [73]. 
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= v h - 








- (xt (1 









for any tensor quantity. Following [74], one can find the scalar gauge invariant variables at first order 
given by 

^=^ + ^{sf 1)a )' , (4.14) 

^e^ + IvV 1 '-^;,, (4.15) 
= fHl) + ( m )' S{[ v (4.16) 
A«=P (1) + (P (0) )'4 ) , (4.17) 

I ( nm (x ll(1, )'A 
with 5^ = I cj i|( - > — - — g-^- I ■ The vector modes are 

(4.18) 

(4.19) 
(4.20) 

Another gauge invariant variables are the 3-current, the charge density and the electric and magnetic 
fields, because they vanish in the background. The tensor quantities are also gauge invariant because 
they are null in the background (see eq.(2.15)). 

4.1 The Ohm's law and the energy momentum tensor 

Using eq. (3.7) the Ohm's law at first order is 

jfWif). (4.21) 

At first order we suppose the conductivity of the medium finite (real MHD) therefore the electric field 
and the 3-current are nonzero. Now the electromagnetic energy momentum tensor is 

(4.22) 
(4.23) 

(4.24) 

where , = | (S 2 + i? 2 ) 5\ — B iJ3 l — E\E l is the anisotropic stresses which is gauge invariant by 

definition cq.(2.15), this term is important to constrains the total magnetic energy because it's source 
of gravitational waves [33]. We see that electromagnetic energy density appear like a cuadratic term 
in the energy momentum tensor, this means that electromagnetic field should be regarded a \ order 
perturbation. 8 Using eq. (3.5) and considering the fluctuations of the matter fields eqs. (4.1) and 
(4.4), the energy momentum tensor for the fluid is given by 

r (%o = -A (1) + (/i ( o))'4)> ( 4 - 25 ) 
T\ SV) o = (/io + Po) fa - 0{ x) - , (4-26) 

7 , (% i = -(Mo + Po)vf ) , (4.27) 



rpO 

1 (em) 




E fi)) , 




rpi 

1 (em)0 


= o, 


(em)i 


o, 




rpi 

1 (em) I 


1 

47T 


[I K, 


+ £(i)) 6\ - 


, TTi(l) 

r 11 ;(em) 



-(fl) 

1 (fi)j 



(A«-(P (0) ) / 5ll ) )^+n'W ) , (4.28) 



7 With 9 ! x!j )T = 0, Xi k)l = and D u = d t d 3 - §5yV 2 . 

8 Thercfore the magnetic field should be splitting as B i = (-B^ + + _B* 2) + .... ), sec [75, 76] 
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where n'^jys the anisotropic stress tensor [77]. The above equations are written in terms of gauge 
invariant variables plus terms as sl^ that depend of the gauge choice. 

4.2 The conservation equations 

The total energy momentum conservation equation Tg. a = can be split in each component which is 
not conserved independently 

n <a = TfJ+Tf^=Q, (4.29) 

where 

Tpg ' M) = Fpaj a . (4.30) 
Using the eqs. (4.25) and (4.28), the continuity equation T$. a = is given by 

(AW)' + 3H (a£> + A«) - 3 (P + mo) + (Po + Mo) VV 1 ) - 3tf (P (0) + M(o))'5^ 

- ((^(o))'5jl ))' + (Po + Mo) (-^V 2 X « + 3P5JI ))' - (Po + mo) V 2 Qx" (1) )' 



i) 
0. 
(4.31) 

The Navier- Stokes equation T" a = is 



V 1 > (M(o)+-P(o)) 1 ' 0(0) +JW aV ( ) J 



32) 



The above eqs. (4.31) and (4.32) are written is terms of gauge invariant variables in according 
to [68, 77, 78], we see that doesn't exist contribution of electromagnetic terms to the conservation 
equations. In [29, 79] the energy-momentum tensor of each component is not conserved independently 
and its divergence has a source term that takes into account the energy and momentum transfer among 
the components of the photon, electron, proton and electromagnetic field = Qp- 

5 Maxwell's equations and the cosmological dynamo equation 

The Maxwell's equations are written as 

V a P^=/, V [7 P Q/3] =0. (5.1) 

Using cq. (5.1) and the pertubation equations for the metric and electromagnetic fields, we find the 
non- homogeneous Maxwell's equations 

^£^=00(1), (5.2) 

e^diBfa = (E l {1) )' + 2HE\ l) + aJ^, (5.3) 

and the homogeneous Maxwell's equation 

B' H1) +2HB { k 1) +e l ld l E^=0, (5.4) 
d l B\ 1] + ^Bf ] + d k B^ = 0, (5.5) 

which is in according to [80]. Now using the above equations eqs. (5. 2), (5. 4), (5. 3) and (5.5) with the 
ohm's law eq. (4.21), we get an equation which gives us the evolution of magnetic field at first order, 
it's called dynamo equation: 



(B^)' + 2HB^+rj 



V x ( V x B {1) - ( P (1) ) - 2PP (1 



= 0, (5.6) 

k 
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with 77 = the diffusion coefficient [4]. The eq. (5.6) is similar to dynamo equation in MHD but 
it's in the cosmological context [3, 4]. This equation has one term that depends on rj diffusion term 
which takes into account the dissipation phenomena of the magnetic field (the electric field in this 
term in general is dropped if we neglect the displacement current). We should have into account that 
77 is a expansion parameter (due to ^ is small), from eq. (5.6) we see that for a finite rj the diffusion 
term should not be neglected. Is important to remark one should take care with the assumption rj = 
to any scale, because it could break to some scale [70]. 

6 Generalization at second order 

Following to [62] is introduced the variable 6^ 2 'T defined by 

4 2) t ee sfr - 2L X (#>r) + L 2 x r Q ; (6.1) 

inspecting the gauge transformation eq. (2.16) one can see that 5^ 2 'T is transformed as 

S§ ) T-6<g ) T = L (T r , (6.2) 

with a = £2 + [£ii^"] an d X is the gauge dependence part in linear order perturbation. The gauge 
transformation rule eq. (6.2), is identical to the gauge transformation at linear order eq. (2.15). This 
property is general not only happens in the case of cosmological perturbation around FLRW metric, 
it's the key to extend this theory to second order 

L [8 2 T] = S [ST, ST] . 

We shall see that first and second order equations are similar however, the last have as sources couples 
between linear perturbations variables. Using cqs. (2.16) and (6.2) we can arrive to the gauge invariant 
quantities at second order. The scalar gauge invariant at second order are given by 



* (2) ee i>W + i (s{l )0 )' + T\0W), (6.3) 
$^ ee 0( 2 ) + iv 2 X (2) - HSf 2) + T 2 (0 (2) ), (6.4) 

A^^^ + ^o^'^+r 3 ^ 2 '), (6.5) 
A( 2 W 2 ) + T 4 (0 (2) ), (6.6) 
Af^Bf 2) +T 4 (0^), (6.7) 
A^^Ef^+T 6 ^), (6.8) 

^ 2 ^«( 2 )+Qx ll(2) )' + r 7 (0 (2) ), (6.9) 
with Sjj^ ee ^JK 2 ) - ( X ",P ^ + T 8 (C (2) ). Vector modes found arc as follows 

^ 2 )=^ 2 ) + (xl (2) )' + r 9 (0^), (6.10) 

^ 2) ^^ (2) -( Xl ±(2) )' + r 10 (O (2) ), (6.11) 

V ( l 2) ^^ 2) +^ 2) +r n (0( 2 )), (6.12) 

B U = B U + - J^ s Mih> ( 6 - 13 ) 

4) = E h) + " j^cVuu- ^ 

j^^,r {2) + r 12 (o^), (6.15) 

n,(f T ee ng>* + ng>» m + t 13 (o^), (6.16) 
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where we calculate the gauge invariant quantities eqs. (6.13) and (6.14) for electromagnetic fields 
because they are gauge dependent at second order. In the eqs. (6.3) to (6.16) we have written the 
gauge invariant quatities at second order. We can see that these variables are similar to these ones 
at first orden but in this case appear sources as T k (JD^) that depend of the gauge choice and the 
couples of terms at first order. The explict calcules of T k {0^ 2 ') can be seen in [61, 74]. 

6.1 The Ohm's law and the energy momentum tensor 

Using eqs. (3.7), (4.5) and (4.6), we get the Ohm's law at second order 



2(7 



s 2 (o { 



2) 



(6.17) 



In this case we see that 3-current has a type Lorentz term and shows couples between first order terms 
that can affect the evolution of the current. Hereafter the functions S^O^) with n G Z gives us the 
gauge dependence. Doing the expansion at second order in the fluid energy momentum tensor, one 
can find the following expression 
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which is in according to [78]. Now consider eq. (5.4) the electromagnetic momentum tensor at second 
order is 
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Using eq. (4.29) the continuity equation is given by 

(AW)' + 3H (a£> + Af) - 3 (V 2 >)' (Po + Mo) 
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and the Navier-stokes equation 
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which is in according to [78, 81]. In this case we see that electromagnetic fields affect the evolution 

(2) (2) 

of matter energy density AJj and the peculiar velocity V\ therefore, these fields influence the 
large structure formation and can leave imprints on the temperature anisotropy pattern of the CMB 
[68, 82-84]. 

7 The Maxwell's equations and the cosmological dynamo at second order 

Using the eq. (4.21), the non homogeneous Maxwell's equations are 

diSL = -iEf^di - 3*W) - Su (O^) + aAf\ (7.1) 



(2) - -<±-^(i)<A I * ' - ov> j - a u iu' j -r u.^ e 

V x B {2) ) 1 = (£ l {2) ) ' + 2H£\ 2) + 2E\ 1) ( 2 ( ' - 6 ( ' 



+ 2 (2*« - 6$«J (V x + a^ 2) + S J ig fO^J . (7.2) 

While the homogeneous Maxwell's equations 

+ djBf + d k B {2) = ~S 16 (o (2) ) , (7.3) 

(Sf )' + 2H (Sf ) + (V x E m ) k = -S? (o^) ; (7.4) 

Where the S k terms carry out the gauge dependence. Using the above Maxwell's equations together 
with the Ohm's law at second order and following the same methodology for the first order case, we 
get the cosmological dynamo equation which describes the evolution of the magnetic field at second 
order 
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+ x 

We find that perturbations in the space-time play an important role in the evolution of primordial 
magnetic fields. In this case we see from eqs. (5.6) and (7.5) dependent on geometrical quantities 
(perturbation in the gravitational potential, curvature, velocity ...). We know that these quantities 
evolve according to the Einstein field equations (the Einstein field equation to second order are given 
in [59, 81]). In this way, Eq. (7.5) tell us how the magnetic field evolve according to scale of the 
perturbation. In subhorizon scale, the contrast density and the geometrical quantities grow, hence, 
the dynamo term should amplify the magnetic field. As a final comment we point out that in order to 
solve the dynamo like equation for the magnetic field is necessary to solve the Einstein field equations 
to the second order together with the conservation equations. 



8 Weakly magnetized FLRW-background 

In this section we consider the background space-time as a single barotropic fluid permeated by a 
weak magnetic field, therefore the energy momentum tensor is formed by two components, the fluid 
and electromagnetic energy momentum tensor [66, 67]. We allow the presence of a weak magnetic 
field, this field has the property B 2 ^ <C /i( ) and must to be sufficiently random, to satisface (Bj) = 

and /b 2 q] \ = (b^'B^-.) =/= to ensure that symmetries and the evolution of the background remain 
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unaffected. 9 Again we work under MHD approximation, thus in large scales the plasma is globally 
neutral, charge density is neglected and the electric field with the current should be zero, thus the only 
zero order magnetic variable is B^ Q s [67] . Bianchi models are often used to describe the presence of a 
magnetic field in the universe due to anisotropic properties of this metric. However, as we are dealing 
with weak magnetic fields, is worth to assume the presence of a magnetic field in a FLRW metric as 
background. 10 Indeed, the authors in [85] found that, though there is a profound distinction between 
the Bianchi I equations and the FLRW approximation, at the weak field limit, these differences are 
reduced dramatically, and therefore the linearised Bianchi equations are the same with the FLRW 
ones. Under these conditions, we find that to zero order the electromagnetic energy momentum tensor 
in the background is given by: 
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The magnetic anisotropic stress is treated as a first-order perturbation due to stochastic properties 
of the field, therefore it does not contribute to the above equations. 11 We can see in eqs. (3.5) and 
(8.1)-(8.3), that fluid and electromagnetic energy-momentum tensor are diagonal tensors, that is, are 
consistent with the condition of an isotropic and homogeneous background [67]. If we consider the 
average magnetic density of the background different to zero, the pcrturbativc expansion at k— th 
order of the magnetic density is given by 
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where at first order we get a gauge invariant term which describes the magnetic energy density 
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one can find that average density of the background field decays as B? Q s ~ fl4 ^ 



At first orden 



we work with finite conductivity (real MHD), in this case the electric field and the current becomes 
nonzero, therefore using the eq. (3.7) and assuming the ohmic current is not neglected, we find the 
Ohm's law 
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The electromagnetic energy momentum tensor is given by 
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where 
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(8.11) 



9 we assume that at zero order the magnetic field has been generated by some random process which is statistically 
homogeneous so that just depend on the time, {..) denotes the expectation value. 

10 In fact, this is what we do every time we study the radiation era. The stress energy tensor of a random radiation 
field is nothing else but the electromagnetic energy momentum tensor without it's off-diagonal terms. 

11 Again, we omit the brackets because the density evolves exactly like field itself [67, 73]. 
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and n' ; ^ m j = i (^Amag + E 2 ^j S\ — BiB 1 — E \E l is the anisotropic stress which appear as a pertur- 
bation of the background, this term is important to constrains the total magnetic energy because 
it's source of gravitational waves [33]. The above equations are written in terms of gauge invari- 
ant variables plus terms as Sn\ which are gauge dependent. Now, using the above equations eqs. 
(5. 2), (5. 4), (5. 3) and (5.5) with the ohm's law eq. (8.6), we arrive to the dynamo equation which give 
us the evolution of magnetic field to first order 
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+ (Vx (B {0) x V (1) )) fc = 0. (8.12) 



When we suppose a weak magnetic field on the background, in dynamo equation appears a new term 
called dynamo term which can amplify the magnetic field, this term depends of the evolution of Vm, 
see eq. (4.32), and also from eq. (4.32), seems likely when matter and velocity perturbation grow the 
dynamo term amplifies the magnetic field, it is a difference with the first approach where the dynamo 
term just appears at second order. Using eq. (3.7) we find to second order 
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In this case we see that 3-current has new terms depending on background magnetic field. Doing the 
expansion to second order in the electromagnetic energy momentum tensor, one can find the following 
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Using the Maxwell equations together with the Ohnrs law to second order eq. (8.13) and following the 
same methodology for the first order case, we obtain the follow relation for the evolution of magnetic 
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fields at second order: 
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(8.19) 



New terms arise because the magnetic density is nonzero on the background, giving origin to terms as 
(V(2) x i?( )) and couples with geometrical terms (2$W + 'I'' 1 ') (V(i) x which could contribute 

to the amplification of the field at small scales. In this approach does not appear fractional orders in 
the variables making a more consistent treatment of the perturbation theory introduced in section 2. 



9 Discussion 



A problem in modern cosmology is to explain the origin of cosmic magnetic fields. The origin of 
these fields is still in debate but they must affect the formation of large scale structure and the 
anisotropics in the cosmic microwave background radiation (CMB) [87-89]. We can see this effect in 
eq. (6.26) where the evolution of ; depends on the magnetic field. In this paper we show that the 
perturbed metric plays an important role in the global evolution of magnetic fields [90]. From our 
analysis, we wrote a dynamo like equation for cosmic magnetic fields to second order in perturbation 
theory in a gauge invariant form. We get the dynamo equation from two approaches. First, using 
the FLRW as a background space-time and the magnetic fields as a perturbation. The results are 
eqs. (5.6) and (7.5) to second order. The second approach a weakly magnetic field was introduced 
in the background space-time and due to it's statistical properties which allow us to write down 
the evolution of magnetic field eqs. (8.12) and (8.19) and fluid variables in according to [67]. We 
observe that essentially, the functional form are the same in the two approaches, the coupling between 
geometrical perturbations and fields variables appear as sources in the magnetic field evolution giving 
a new possibility to explain the amplification of primordial cosmic magnetic fields. One important 
distinction between both approximations is the fractional order in the fields which appears when we 
consider the magnetic variables as perturbations on the background at difference when the fields are 
from the beginning on the background (section 8). Although the first alternative is often used in 
studies of GWs production in the early universe [75, 76], the physical explanation of these fractional 
orders is sometimes confused, while if we consider an universe permeated with a magnetic density 
from the background, the perturbative analysis is more straightforward. Further studies as anisotropic 
(Bianchi I) and inhomogeneous (LTB) models should be addressed to see the implications from the 
metric behaviour in the evolution of magnetic field and relax the assumption in the weakness of the 
field. 
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